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Abstract 



A striking contradiction between lattice short range static potential (n/ = 0) and 
standard perturbative potential, observed in Ref. [18], is investigated in the framework 
' of the background perturbation theory. With the use of the background coupling a b (r) 

which contains the only background parameter - mass m b , fixed by fine structure fit in 
bottomonium, the lattice data are nicely explained without introduction of exotic short 
in ' range linear potential with large " string tension" cr* ~ 1 GeV'^ . A significant difference 

between asir) and standard perturbative coupling av{r) is found in the range 0.05 fm 
Q>^ , ^ ^ 0.15 fm, while at larger distances, r > 0.3 fm aB(r) fast approaches the freezing 

O I value aBioo). Some problems concerning the strong coupling properties at short and 

^ ■ long distances are discussed and solutions are suggested. 

o , 

■ 1 Introduction 

Q^l The property of freezing of the strong coupling constant avi''^) at long distances is widely used 
^ ' in QCD phenomenology On the fundamental level this phenomenon has been studied 

in two different theoretical approaches 10]- 0- Iii the case of the static potential the freezing 
of the coupling avir) suggests that ai/(r) is approaching a constant = av{r oo) at 
rS I relatively long distances while at small r it manifests the property of asymptotic freedom. 
c3 I Both characteristic features of the static potential were widely used in hadron spectroscopy. 

However, it was realized that the asymptotic freedom behavior does not practically affect 
hadron spectra being important mostly for a wave function at the origin. On the contrary, 
the choice of (5y(r) as a constant at all distances, i.e. (5y(r) = a, appears to be a reasonable 
approximation and gives rise to a good description of meson spectra both for heavy quarkonia 
[H, |ll| and for heavy-light mesons |T^. Also in lattice QCD this choice gives a good fit 



to lattice static potential at the distances above 0.2 fm. 

Therefore the question arises why this simple approximation, ay(r) ^ a, works so well 
even in the case of bottomonium where the sizes of low-lying levels are not large, the charac- 
teristic radius Rch ^0.5 fm. To answer this question one needs to clarify another problem. 
Namely, to find out the precise freezing value of the vector constant in momentum and co- 
ordinate spaces, and to define the distances r where the difference between ay(r) and is 
becoming inessential and therefore the approximation C(v{r) = a {a ^ Ofr in general case) 
gives a good description of hadron spectra and other physical characteristics. 
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This problem will be discussed in the present paper in the framework of background 
perturbation theory (BPTh) and it will be shown that the background coupling asif') ap- 
proaches its freezing value already at rather small distances r > 0.4 fm. Here it is worthwhile 
to remind that in most calculations in coordinate space the " average" value a is usually taken 
in the range: 0.35 < a < 0.45 while in momentum space larger critical values were used 
(acr = avil = 0)), for example, in Ref. [Q acr = 0.60 and for the Richardson potential 



On 



2n 



0.7 (n/ = 3) 0, 13|, i.e. the difference between ~ 0.6 -t- 0.7 in momentum 



space and a ~ 0.40 ± 0.05 in coordinate space is essential. However, this difference was not 
confirmed by the analysis of the QCD coupling in background fields where the vector 



coupling constants asiq) and aB{r) were found to have the same asymptotic value: 

= 0) = aB(r oo) = afr. (1) 

This equality takes place also for the phenomenological coupling taken as a sum of Gaussians 
in Ref. §. 

In lattice measurements of the static potential at long distances the freezing phenomenon 
is also seen, however, existing lattice data have not clarified our knowledge about a^. As 
shown in Refs. ||T5[, |T^ lattice static potential at r > 0.2 fm can be parameterized with 



a good accuracy by Cornell potential with rather small a (in lattice notation |a = e). In 
quenched approximation (nj = 0) the fitted lattice values of a ~ 0.20 -hO. 24 (e = 0.27-i-0.32) 
turned out to be small so that in some cases there appears a discontinuity of the vector 
coupling constant at the matching point, rmat ~ 0.2 fm |^. But if dynamical fermions are 
introduced, in lattice QCD the fitted value of a {rif = 2,3) was found to become larger: 
a = 0.30 (e = 0.40) [jl^ still being less than in phenomenological models. 

Another problem concerns behavior of av{f) at short distances. The most interesting 
and unexpected results were obtained in lattice measurements of the static potential at short 
distances, 0.05 fm< r < 0.15 fm |]18|, where a large difference between lattice and two-loop 
(one-loop) perturbative potential was observed, yielding discrepancy about 100% at the point 
r = 0.15 fm. In Ref. this large difference was parametrized by a short distance linear 
term a*r with the slope a* ~ 1 GeV^. 

This effect will be explained in our paper. To this end the strong coupling constant in 
background fields asif) will be calculated and the influence of the background mass ttlb will 
be shown to become essential already at rather short distances. We have found that there is 
no need to introduce an additional exotic linear potential a*r as in Ref. [18]. Resulting static 
potential VB(r) appears to have an effectively linear term in a good agreement with lattice 
data. In our calculations no fitting parameters are introduced: the value of the background 



mass niB = 1.0 GeV is taken from fine structure analysis in bottomonium |TT[ while the 
QCD constant Aj^g^Uf = 0) is considered as a well established number and taken from 
Ref. PI . 



2 The strong coupling constant aB{q) in background 
field theory 

The perturbative static potential is used to define a coupling constant av{q) in \^-scheme: 

VA,) = -4vrC.^, (2) 
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where = q^. Recently the renormahzed av{q) was calculated in two- loop approximation 
[ PD| , |2T[]. In coordinate space the static potential can be defined as the Fourier transform of 
Vpiq), 

Vp{r) = J ^\/p(g)exp(zqr), (3) 
which gives rise to the simple relation between the coupling constants: 

~ r \ 2 /-^ singr 

(^v{r) = - / dq av{q), (4) 

TT JO q 

if the following definition for the coupling in coordinate space ay(r) is used, 

Vp{r) = -Cf — - — • (5) 

However, the Fourier transform of the perturbative coupling av{q), Eq.(|D, does not exist 
in a strict sence because of Landau pole singularity. To escape IR divergency the expansion 
of av{q) in the perturbative series at large is usually made |2y, |21|, but the resulting 
expansion is valid only at short distances. 

Here we suggest to obtain the static potential in coordinate space with the use of the 
coupling in BPTh, where the vector coupling constant apiq) in momentum space is defined 
at all momenta and has no singularity for g^ > |0] . Then the potential in momentum space 
can be written as in Eq.(0), 

Vsiq) = -CpArr^^. (6) 

In Eq. (6) and below we consider the infiuence of background vacuum fields only on 
Coulomb-type interaction. In presence of background fields the QCD coupling is modified 
so that it depends on the combination (g^ -|- m^) instead of as it is in standard pertur- 
bative theory The mass ms is a background mass which is characteristic for a process 
considered. In two-loop approximation the running background coupling is 

= (7) 



where one-loop expression is given by 



a4"(,) = (8) 

with 

*B=ln^. (9) 

The condition ttib > Ay is assumed to be satisfied under the logarithm (P) to guarantee the 
absence of Landau pole; this condition is always valid for the numbers Ay and rriB used in 
our calculations (see the numbers in Eq.(|25D). 
In Eq.(0) 

/?o = ii-^%;A = i02-i;. (10) 
First we discuss the most important properties of the background coupling ^^(g). 
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The background mass is not an arbitrary parameter. It can be calculated in the 
framework of BPTh or in lattice QCD. As was shown in Ref. |^ the background mass 
in the case of the static potential is defined by the difference of two-gluon and one-gluon 
hybrid excitations and can be extracted from the corresponding level differences of hybrids 
ccg,bbg. In Ref. this mass rriB was evaluated to be 1.0^1.2 GeV. For other processes, 
like e'^e^ —>■ hadrons, in general case the background mass rriB may be different p|. The 
appearence of the mass in Eq.(9) is similar to the case of QED where a has the mass of 
e^e~ pair under logarithm. 

It is of interest to notice that the analytical form of aB{q) (0) coincides with that obtained 
in a picture when gluon is supposed to have an effective mass rrig inside gluon loop. Therefore 
in Refs. Cisil) was taken as a function of (g^ + 4mp, i.e. the double effective gluon mass 



2mg plays a role of the background mass rriB (see a discussion in Ref. 0]). 

In our calculations here the value of ttlb will be fixed at rriB = 1.0 GeV taken from the 
fit to fine structure splittings in bottomonium [pT] . 

At large momenta, ^ m%, the background coupling goes over into the standard per- 
turbative expression av{q)- Therefore the QCD constants in PQCD and A^ in BPTh 
must coincide, in any case it is true for the number of flavors n/ = 5. As can be directly 
calculated with the use of matching procedure they are also equal for Uf = 4. 

The QCD constant Ay, entering the coupling av{q) in Eq. @, can be expressed through 
Aj^(n/) in MS renormalization scheme p4| : 



with 



111 



31 10 



(12) 



At present the values of A^^"* are well established for the number of flavors nj = 5: A^^ 



208±|| MeV p5| and also for Uf = due to analysis in lattice finite size technique [ITS 



MS 



A. 



(0^ 
MS 



602(48) 



(13) 



where tq denotes the Sommer scale. With the use of Tq = 2.5 GeV , taken in most lattice 



calculations 17, 18 , one obtains 



A. 



(01 

MS 



241 ± 19MeV, 



then from Eq. (|TT]) the QCD constant in the V^-scheme is 



(14) 



A 



(0) 



385 ± 30MeV. 



In three-loop approximation the background coupling 



«g^(?) =<(?)<! 1- 



+ 



Ib PotB 

contains the term including coefficient, which depends on a renormalization scheme and 
was calculated in Refs. ||20|, ^ (the coefficients Po,Pi do not depend on the renormalization 



{IntBY-lntB - 1 + 



Pi 



(15) 



(16) 
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scheme); this coefficient 

/?2'' = /?r-«lA + («2-«?)/30. (17) 

Here ai is defined by Eq.(|l2|) and 

2857 5033 325 , 
^2 2 18 ^ 54 ^' ^ ^ 

a2 = 9 + 47r2 + —(3 + nl- 

V 162 4 3 7 81 ^ 

3 /1798 56 . _A 2 /55 



In Eq.(^9|) ({3) = 1.202057 denotes the Riemann (^-function. In quenched approximation 

a2(n/ = 0) = 456, 7488, p^{nf = 0) = 4224, 1817, (20) 

i.e. coefficient turns out to be about 3 times larger than 13^^ {nj = 0) = ^y^. As a result 

the third order correction in a^sio) much larger than a^'*(g) and a^sio)- 

It is easy to find the first order correction to the perturbative coupling ay{q) which comes 
from the expansion of ^^(g) in powers oim\/q^. In two-loop approximation 

«il(large q) = 4\q) - a\^\q)^^; (21) 



appr \ 

q' -L 



with 

K y Poy Ay 

This approximation appears to be valid only at g > 2 GeV (with the accuracy < 10%). 

The behavior of asiq) in IR region. The freezing value of a;^''(g) can be easily obtained 
from Eqs. (0)-(^, in particular, for two- loop coupling 

(2) (2). 2 47r J /?i Intol .^on 

Po^o t Po f-o J 

with 

«L'^ = 7T-; to = tB{q' = 0)=\n-^. (24) 

In what follows the notation a^'' = ag''(g^ = 0), as in potential models [^, is also used. 
The parameters and Ay, present in ts, are considered to be fixed: the value of is 
taken from fine structure analysis of 2P- and IP- states in bottomonium while the value 
Ay\nf = 0) is taken from lattice data and given by Eq.(|T5|), 

mB = 1.0GeV, AS!?^(n/ = 0) = 385 MeV. (25) 

We suppose here that for Uf = 3 the constants Ay"* and Ay"* are approximately equal, 

Ay(nf = 3) ~ Ay(n/ = 0) = 385 ± 30 MeV (26) 
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Then the following critical values can be obtained from Eqs. (|7|) , (|T5|) : 

= 0.598; a^?) ^ o.428; ag^ = 0.805 (n/ = 0) (27) 

ail'^ = 0.731; a^'^ = 0.536; = 0.972 (nj = 3) (28) 

As one can see from Eqs. (|27D , (|28|) the third order coupling turns out to be about 90% 
(uf = 0) and 80% (n/ = 3) larger than a^\nf) because of large P2 coefficient (20). Such 
a value of a^\nf) appears to be too large and not compatible with the effective Coulomb 
constant of the static potential, ^ 0.5, needed to describe heavy quarkonia [|l], |^ and 
heavy-light meson spectra [|12|. Such large value of Coulomb constant is also not observed 
in lattice calculations of the static potential at large distances [15,16]. 

Presumably it means that since the perturbative series is an asymptotic one, it should 
be truncated after the second-loop term. 

Recently an "analytical" perturbation theory was elaborated in papers [9]. It was shown 
there that the modification of singularities of perturbative coupling by power terms allows 
to work accurately in two-loop approximation. The third order contribution becomes nu- 
merically inessential and e.g. for e"'"e~-annihilation is about 0.5%. 

Therefore in what follows the third order term will be omitted and a posteriori our 
phenomenological analysis with the use of background coupling will demonstrate that the 
two-loop approximation is sufficient to describe the lattice data at small distances. 

The background coupling a^'' (q) in two-loop approximation turns out to be rather close 
to the phenomenological aph(g) which is successfully used in hadron spectroscopy. For com- 
parison apii{q) will be taken from the well-known paper of Godfrey and Isgur 

aci = 0.25 exp(-g2) + 0.15 exp(-0.1g2) + 0.20 exp(-0.001g2) (29) 

with q in GeV and acr = 0.60. 

In Fig.l this phenomenological coupling is compared to the background coupling (g) 
for n/ = 3. Here, as in Eq . (^6]) , it is supposed that Ay(n/ = 3) = Ay(r2/ = 0) and for 
Av{nf = 3) two values, are taken: 

(A) Avirif = 3) = 385 MeV 

(B) Ay(n/ = 3) = 410 MeV. (30) 

These values of Ay-* do not contradict those which are commonly used in MS renormalization 
scheme, and give rise to as{Mz) = 0.118 ± 0.001. The connection between A^^ and A^y^^ 
is given by Eq.(ll). 

In Fig.l the couplings aciiq) (solid line), a^Ho') i'^ two-loop approximation with Ay = 

(2) 

385 MeV (dashed line), and as (q) with Ay = 410 MeV (dash-dotted line) and = 
1.0 GeV are shown. At the momentum q = rhc = 1.3 GeV [rhc is the running mass of 
c quark) the matching of the couplings was done with the following result for the QCD 
constant: Ay(n/ = 4) = 0.325 MeV in the case A and Ay(n/ = 4) = 351 MeV in the 
case B. As is seen from Fig.l the background coupling ^^''(Q') iii the case B (and to some 
extent in the case A) appears to be very close to the phenomenological coupling q;gi(?); the 
difference between them is less than 5% at small q < 1.3 GeV and less than 2% in the range 
1.3 < g < 4 GeV. So one can expect that with the use of the background coupling as good 
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description of low-energy experimental data can be obtained as in Ref. with the use of 
the phenomenological coupling. From here our estimate of the freezing value is about 0.53-i- 
0.60, and it is interesting to look at the expansion of the background coupling asiq) near 
the freezing point g = 0: 

«ippr (small q) = a^^ - a^^ — ^ , (31 ) 

In -j^ 



where a^^^^^cr'' ^^e fixed numbers defined by and Ay (see Eqs. (PB|), ([2^)). This 
approximation appears to be valid only in very narrow range of small momenta q,0 < q < 0.4 
GeV, where the difference between a;^^(g) and a^p^pj. (small q) is less than 5%; but it already 
reaches 22% at g = 1.0 GeV. 

Heavy-quark initiated jets can be successfully described at small momenta if the following 
assumption is made about an effective coupling constant a;efr(g) in infrared region 0: 

/■2 GeV fy „(fj] 

J2(fit) = (2 GeV)-^ / dq^^ = 0.18 ± O.Ol(exp) ± 0.02(th). (32) 

JO 71 

As was shown in Ref. this number does not depend on the form of aes{q) assumed in the 
fit. Our calculations of the integral (|32|) with the background coupling a^B\q) from Eq.(|^) 
with Uf = 3, Ay = 410 MeV, tub = 1.0 GeV give 



r2 GeV (^Q , , 

J^^as) = (2 GeV)"^ / — «b (?) = 0-134 (33) 

JO TC 



dq 



i.e. this number is by 26% smaller than J2(fit) in Eg . (p2|) . The same integral calculated with 
the phenomenological constant aQi{q) (^) is also by 22% smaller than J2(fit) (the central 
value): 

J2(aGi) = 0.14 (34) 
and very close to our number (P^l). One should notice here that the large number (^) in 
Ref.[6] could be connected with the large fitted value of a^^(M^) = 0.125 ± 0.003(exp) = 
0.004(th), used in their paper, while now the average a^^{Mz) = 0.118 ± 0.001 is accepted 

In conclusion we give our predictions about the freezing values in momentum space: 





= 0.598, 




= 0.428 


{nf = 0, 


a(0) 


= 385 MeV), 


a« 


= 0.731, 




= 0.536 


{uf = 3, 


a(3) 
Ay 


= 385 MeV), 




= 0.783, 




= 0.582 


{nj = 3, 


a(0) 
Ay 


= 410 MeV) 



(35) 



3 The background coupling ^^(r) in coordinate space 

From the explicit expressions of ■* (g) it is evident that the coupling a^^^ (g) is well defined 
at all momenta g^ > if the condition > Ay is satisfied. Therefore the Fourier transform 
can be used to define the static potential in coordinate space over all distances: 

Vb r = -Cf = -Cp^'n / — —e ^ ——^ 36 

r J q^ Utt r 
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From here the relation similar to Eq.(H) follows: 



^ , , 2 f°° , sin qr , , 2 sinx , , , 
dByf) = — cLQ ctByq) = — dx aB{x/r), [Si) 



TT Jo q n Jo X 

where now the background coupling asix) depends on the variable 



tB{x)=ln J (38) 

IVyf 



This integral (^) cannot be taken analytically even in one-loop approximation and was 
calculated numerically in n-loop approximations {n = 1,2) with the use of the parameters 
(p5|). The behavior of background coupling a^^\r) and a^^\r) is shown in Fig. 2 in the range 
< r < 1.4 fm. 

From Fig. 2 one can see that two-loop background coupling in coordinate space is ap- 
proaching the freezing value at relatively short distances, r > 0.4 fm, and the values of 
a^^\r) = 1, 2) at the Sommer scale tq ^ 0.5 fm are following, 

ag^(ro = 0.5 fm) = 0, 574; ag^(ro = 0.5 fm) = 0.404 (39) 

It is of interest to notice that the two-loop coupling at the distance tq practically coincides 
with the number a = 0.39 widely used in Cornell potential while the one-loop coupling is 
too large. This fact that two-loop background coupling is almost constant already at r > 0.4 
fm can be considered as an important argument in favour of the choice OBif) =const in low 

— (2) 

energy spectroscopy. In Fig. 3 (r) is shown for two different values of the background 
mass: tjib = 1.0 GeV (solid line) and rriB = 1.1 GeV (dashed line). As seen from Fig. 3 the 
difference between them is becoming essential already at r 0.3 fm, being about 10% over 
all distances r > 0.3 fm; their freezing values are af^\mB = 1.0 GeV)=0.428, af^\mB = 1.1 
GeV)=0.382 (in both cases = 385 MeV). 

The freezing value in coordinate space turns out to be just the same as in momentum 
space, i.e. for Uf = 3 and Uf = [Ay = 385 MeV) they are given by the numbers from Eq. 
(pSl). This property is true also for the phenomenological potential used in Ref.|p|: 



aciiq = 0) = «Gi(r -> oo) = 0.60, (40) 

since in coordinate space the coupling aGi{r) corresponding to q;gi('?) in Eq. (|2^) is 

aGi(r) = 0.25$(2r) + 0.15$(1.581r) + 0.20$(15.811r), (41) 

where ^{z) is the error function. Thus for three flavors the phenomenological value ~ 0.6 
was found to be a bit larger than our number = 0.54 (see Eq. (^) for = 3). 

With a_B(r) calculated above we can compare the background potential VB(r) (37) to the 



lattice static potential from Ref.[T^. Here we are mostly interested in short range potential, 
in particular, in the influence of background mass on its behavior. The properties of 
asif) at small r will be considered in the next Section. 



4 (r) at short distances 

Recently very precise lattice measurements of the static potential at short distances were 
presented [|^. Having these data one has a unique opportunity to compare theoretical 
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predictions about the background coupling and the potential VB^r) with precise lattice data 
in quenched approximation. We remind that our calculations of a^^\r) (n = 1,2) were 
done without any arbitrary parameter: A^^ = 385 MeV {uf = 0) was fixed from lattice 
data (Eq.(|l^)) and = 1.0 GeV from the fine structure sphtting of IP and 2P states in 
bottomonium. In Fig. 4 (5^''(r) is compared to the perturbative running coupling a[^''(r) at 

~ (2) 

r < 0.12 fm, ay {r) is given by the expression (rA/j ^ 1) 



«V ('^J = «y W i 1 - 7^2 r (42) 

Po Vr 



with 



and the following prescription for the value of QCD constant A/j(n/) p^, [p!3[] 

Ai"^)=Al7^)exp7i. (44) 

In Eq. ( ^41) 'Je is the Euler constant (7^; = 0.5772157), and in quenched approximation A^^-* 
is given by the number ([T^), therefore for nf = 

Ag^ = 684 ± 53 MeV (45) 
In our calculations below we take the number 

Ag^ = 686 MeV, (46) 

which corresponds to Ay"* = 385 MeV according to the relation (H). 

~ (2) 

The numerical comparison of the "exact" background coupling (r) (Av = 385 MeV, 
Uf = 0) and the corresponding perturbative ay\r) with A^'* from Eq. (|46|) is presented in 
Table for the distances in the range 

0.002 fm < r < 0.15 fm. (47) 



Table The background coupling (r) (Ay = 385 MeV, 
niB = 1.0 GeV) compared to perturbative coupling ay\r) with A^"* = 686 MeV. 



r in fm 


4'^(r) 


«l/'^(r) 


r in fm 


4'^(0 




0.002 


0.0964 


0.0924 


0.024 


0.1757 


0.1667 


0.004 


0.11095 


0.10505 


0.030 


0.1880 


0.1807 


0.006 


0.1216 


0.1209 


0.0355 


0.1988 


0.1943 


0.008 


0.1303 


0.1221 


0.041 


0.2085 


0.2079 


0.012 


0.1446 


0.1352 


0.049 


0.2204 


0.2264 


0.016 


0.1564 


0.1465 


0.057 


0.2311 


0.2455 


0.020 


0.1666 


0.1569 


0.063 


0.2384 


0.2605 
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One can see that the difference between these two couphngs, 

Aag)(r) = 5g)(r)-4^)(r), (n = l,2), (48) 

has several prominent features. 

First, at very short distances, r < 0.04 fm, the correction Aa^ (r) is positive, i.e. 

agV) > c^Pir) (r < 0.04 fm), (49) 

and relatively small. It is about 6% at r = 0.02 fm and still remains ~ 4% at much smaller 
r = 0.002 fm so that aP{r) approaches the perturbative running coupling ay\r) rather 
slowly. In one-loop approximation this correction was calculated analytically in Ref.[14]: 



^""''•'- 6ft|-HA.r)|3 - IKAvr)|»l. (50) 

It is positive and less than 5% only at very short distances, r < 0.007 fm. The important 
feature of Aa^^\r) is that it does not depend on the background mass in limit r — > 0. 
Secondly, the values of a^''(r) and ay''(r) coincide at the point r = 0.041 fm, i.e. 

Aag^(r = 0.041 fm) = (51) 

At bigger distances, in particular, in the range 

0.04 fm < r < 0.15 fm (52) 

this correction is negative and fast growing, e.g it is 13% at r = 0.07 fm, already 36% at 
r = 0.10 fm and reaches 100% at r = 0.14 fm although all these points lie rather far from 
the Landau pole, rpoie = 0.29 fm. 

The explanation why the perturbative coupling is essentially larger than the background 
coupling at rather small r was given in Ref . . It was shown there that in coordinate space 
the QCD constant Ar can be defined as a constant (given by Eq.(45)) only at very short 
distances while in the transition region (E^) the role of QCD "constant" plays a function 
Aij(r) dependent on the distance: 

Ar ^ AR{r) ^ Ay exp(7^ + ^"^.f ) = Aiiexp(-mBr + -m|r^) (53) 

k=i '^■'^ ^ 

Then with the use of the function A/j(r) the perturbative coupling ay\r) reproduces aP{r) 
in the range (|52D with the accuracy better than 5%. Actually, this approximation (|53[) can 
be used only at distances 

2 2 

rriBr — < 'Je, or r < 0.15 fm, (54) 



In what follows the region (^J) is called as the transition region. By direct calculations of 
the integral (^) it can be shown that at longer distances A/j(r) is approaching the vector 
constant Ay fl^: 

AR^Ay at r>0.15fm. (55) 
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5 Static interquark potential 



Knowing the differences Aa^''(r) = a^^\r) — ay\r) one can calculate the corresponding 
differences AV^\r) between the background and perturbative static potentials: 

AV<»'M = -i^^^^. (n = 1.2). (56) 

Numerically calculated AVg^'^'*(r) are shown in Figs. 5,6 correspondingly. We observe 
the linear rise of potentials illustrated by the tangents (dashed lines). Tangent slopes cr^''(^) 
are defined as 

dAVlp\r) 

'(r) = . (57) 



We see from the figures that 



*(1) 



1.20 GeV^ 0.05 f m < r < 0.12 fm, (58) 



^ Q_g7 Qgy2^ g 03 < r < 0.09 fm. (59) 

This large difference means that the perturbative avi'f') fails to describe lattice data in this 
region. Nevertheless, we try here to compare our results (|58| ) ,(|59| ) with the lattice ones from 
Ref.[18]. After perturbative potential subtraction from lattice static quenched potential, 
lattice data in the region 0.03 fm < r <0.15 fm were parameterized by linear potential with 
the slope 

al, = (1.20 ± 0.36) GeV^ 0.03 fm < r < 0.15 fm. (60) 

One can see that crjt^^ corresponds to a"^^'^^ within one standard error. We conclude that this 
large linear slope of potential is well explained using ^^(r) instead of av{r). In Fig. 7 we 
compare the lattice static potential from Ref.[18] at the distances 0.05 < r/ro < 0.45 with 
the potential 

V{r) = --^^^^ + ar + C, (61) 
3 r 

where a = 0.2 GeV^ and C = —253 MeV (shown by solid line). This potential includes 
background perturbative potential Vg (r) and linear confining potential or. Constant C 
corresponds to the quark self-energy. We observe that this potential describes all lattice 
data remarkably well. 

At the Fig.7 1-loop and 2-loop perturbative potentials, calculated in Ref.[18], are also 
shown. The difference between them and lattice points is large and explained above. The 
1-loop perturbative + linear potential with the slope a*g_^ is also shown. It describes lattice 
points up to r/ro = 0.35, but fails to describe the rest of data because of the Landau pole of 
ay(r). On the contrary, V{r) (Eq. (|6lD ) not only describes all the presented in Fig.7 lattice 
data in the region 0.05 < r/ro < 0.3, but also all lattice data available up to r = 3ro with a 
reasonable accuracy. 



What is the physics of the linear part of V{r)? In the framework of BPTh it was 
shown that the linear confining potential starts from the quark distances close to the gluonic 
correlation length Tg. From the lattice data Tg ^ 0.2 fm and Tg = 0.12 0.15 fm 



[p9|] . At smaller distances one needs to take into account an interference of perturbative 
and nonperturbative effects As was shown in Ref.||30|, the interference potential at 
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r < Tg behaves like a linear one with the slope about a, while nonperturbative potential 
is proportional to and small in this region p7[. At distances r > Tg, the interference 



interaction vanishes, while the nonperturbative potential becomes the linear one with the 
same slope a. As a result the sum of these potentials may be close to a linear one with 
the standard slope a over all distances. This linear behavior of nonperturbative potential 
at short distances is important for the fine structure fit in charmonium [10]. Note also that 
the linear potential V^^^\r) is consistent with the "short string" potential generated either 
by the topologically defined point-like monopoles or infinitely thin P-vortices within the 
Abelian Higgs model [^. The distinguishing feature of last potential is that it does not 



change its slope at all distances. We leave the detailed numerical analysis of this effect to 
the subsequent paper. 



6 Conclusions 

In our paper the strong coupling a 3(1^) in coordinate space, deduced in BPTh, is investi- 
gated and the corresponding perturbative potential is compared to short range lattice static 
potential from Ref. [18]. The following prominent features of (5b('") are observed: 

1. The background coupling attains the asymptotics of standard perturbative coupling 
Q!y(r) only at very short distances, r < 0.04 fm, where the QCD constant A^^ = Avexp'jE 
and Ay, the QCD constant in the ^-scheme is considered to be a well established number 
in quenched approximation. 

2. At larger r a function Aji{r) plays the role of the QCD constant which at short 
distances, m^r < 'Je, is approximately given by 

2 2 

Afi(r) ^ A^jexp(-mBr + ), (62) 

where uib = 1.0 GeV is the background mass fixed by fit to fine structure splittings of IP 
and 2P levels in bottomonium. In fact the condition 

m^r — -m|r^ < 7^; (63) 

defines the narrow transition region: 0.05 fm < r < 0.15 fm, where ^^{r) decreases almost 
twice, from the value at r ^ 0.05 fm to the number close to Ay at r=0.15 fm. 

3. The static potential Vsir) in BPTh, defined through as(r) in usual way, Eq. (36), is 
in a good agreement with quenched lattice static potential measured at short distances [18]. 

4. The specific behavior of asir) in the transition region produces a linear rise of the 
difference A\^(r) = VB(r) — Vp{r) with the slope a* ~ 1 GeV^ in accord to lattice data. 
Moreover, we have obtained a very good agreement with the lattice data using the sum 
VB{r) + ar with a = 0.2 GeV^. 

5. At distances r > 0.2 fm the function Aji{r) turns out to be almost constant; A/j(r) ^ 
Ay over all distances, so that at r > 0.3 fm the coupling aB{r) fast approaches the freezing 
value as (00). This fact helps to understand why the use of ay(r)=const in static potential 
of quark model [3] appears to be a good approximation in hadron spectroscopy. 

6. The freezing value of the background coupling coincides in momentum and coordinate 
space and this statement can be checked in different processes in low energy QCD. Our 
estimate for afr is 0.53-0.60. 
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Figure 1: The behavior of the background couphng Q;^^(g') for Ay' = 385 MeV (dashed 
hue) and for Ay"* = 410 MeV (dash-dotted hne) compared to the phenomenological couphng 
Q!Gi(?) (sohd hne) taken from Ref.[3]. 
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Figure 2: The 1-loop Q;B^(r) (sohd hne) and the 2-loop Q:^''(r) (dashed hne) background 
couphngs in quenched approximation; in both cases Ay' — 385 MeV, ijib = 1.0 GeV; 2-loop 
asymptotics, a^^' — 0.428, is shown. 
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Figure 3: The background coupling Q!^^(r) in coordinate space for two values of the back- 
ground mass niB- ruB — 1-0 GeV (sohd line) and = 1.1 GeV (dashed line); in both cases 
= 385 MeV. 
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Figure 4: The background coupling a^''(r) with Ay*-* = 385 MeV, ms = 1-0 GeV (solid 
hne) compared to the perturbative ay\r) with Kr = 686 MeV (dashed hne) at short dis- 
tances. 
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Figure 5: The difference AVg^^ between tlie background and perturbative potentials in one- 
loop approximation (solid line). The tangent with the slope a*^^"^ — 1.20 GeV^ is shown by 
the dashed line. The parameteres are the same as in Fig. 4. 
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Figure 6: The difference AVg^'' between the background and perturbative potentials in two- 
loop approximation (sohd hne). The tangent with the slope a*^'^^ — 0.87 GeV^ is shown by 
the dashed line. The parameteres are the same as in Fig. 4. 
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